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I. INTRODUCTION 

Chromatin is an ubiquitous protcin-DNA complex ob- 
served in chromosomes of ah eukaryotic organisms and 
strikingly conserved during evolution [Q |Q It en- 
sures DNA compaction during the mitosis and plays a 
key but still unclear role in all the biological functions 
involving genomic DNA: replication, transcription and 
repair. For instance, the chromatin degree of compaction 
is acknowledged to regulate, somehow, transcriptionally 
active regions ||]. 

The general issue taken up in this paper is to under- 
stand the interplay between the mechanical properties of 
the fiber and its biological functions. We aim at under- 
standing quantitatively the grounds for existence of the 
so universal and so conserved features of the assembly 
of linkers and nucleosomes forming the chromatin fiber. 
We claim that part of the answer lays in the peculiar 
mechanical properties of this assembly. 

The typical length scale of chromatin fiber (diameter 
about 30 nm) is large enough to allow a mechanistic view- 
point: describing the elastic properties of the chromatin 
fiber is nothing but a problem of spring mechanics 
Nevertheless, the architecture of this "spring" is much 
more complex than a simple helical coiling and we expect 
that the detailed structural features of the chromatin as- 
sembly still strongly influence the behavior at the fiber 
scale. We thus investigate the specific elastic behavior 
associated with the chromatin structure and its sensitiv- 
ity to the structural parameters of the fiber, singled out 
within a basic geometric modeling of its assembly. 

We underline that fiber elasticity is involved in two, 
quite different, issues. A first one is to describe the lin- 
ear elastic response of the fiber to global stresses, i.e. a 
force and a torque applied at its ends. This issue refers to 
micromanipulations in which a single chromatin fiber is 
pulled (and possibly will be twisted) 1^ . Our theoretical 
study provides a framework to interpret the experimen- 
tal results, in particular to predict force-extension curves 
in terms of geometric and mechanical parameters of the 
underlying model of fiber. Our approach largely extends 
preliminary results presented quite recently by Schies- 
sel et al. 0, since it gives the exact analytical value of 
the complete set of elastic constants of the fiber. Com- 
parison with the observed curves allows to extract small 
scale information, either about the microscopic structure, 
either about the interactions involved, either about pos- 
sible conformational or structural changes. 

However, except during anaphase, when sister chro- 
matids are separated by the mitotic spindle [||, such 
"macroscopic" stresses are not encountered in vivo at 
the chromatin fiber level; micromanipulations of the fiber 
and associated force measurements may nevertheless un- 
ravel physical parameters involved in biological mecha- 
nisms, such as the elastic constants. An issue directly 
relevant to the in vivo functioning of chromatin is to de- 



scribe the response of the fiber to local, internal stresses 
as those created by intercalators, groove-binding pro- 
teins, or any induced change in the fiber assembly or 
small-scale structure. Our theoretical approach also gives 
a framework to such studies. Indeed, to be solved, both 
issues require to relate the structure and mechanics at the 
DNA scale and those at fiber scale, which is the scope of 
the present work. 

The chromatin scale is precisely the scale of nanome- 
chanics; at this scale, we expect a strong and direct in- 
terplay between the biological functioning, monitored by 
various enzymes, and the mechanical properties of the 
substrate, here the chromatin fiber. We thus believe that 
a mechanistic approach is well-suited to evidence pos- 
sible mechanisms for the fiber decondensation prior to 
transcription, for the connection between enhancer and 
promoter DNA regions during transcription, for the nu- 
cleosomal DNA site exposure allowing protein binding at 
specific sites or for the ejection of nucleosomes presum- 
ably required during replication 

The core of the study is to describe how the linear 
elastic behavior of chromatin fiber originates from the 
elastic properties of linker DNA (i.e. naked DNA con- 
necting the nucleosomes). In Section II, we introduce 
a mechanical model of chromatin fiber, quite similar to 
the so-called two-angle model ||T^; it incorporates micro- 
scopic data such as DNA structure and elastic properties 
lll^ as well as nucleosome cristallographic data jl^] 
into an analytically tractable assembly. Its implemen- 
tation gives the geometric properties of the chromatin 
fiber, presented in Section III with a special mention to 
symmetry properties and to the quantities relevant to the 
mechanics of the fiber. Section IV is devoted to the elas- 
tic properties of the fiber. The first step of their study 
is to relate the distribution of stresses along linker DNA 
to the "macroscopic" stresses (force and torque) applied 
to the fiber ends. This allows to compute analytically 
the elastic constants describing the linear response of 
the fiber from the knowledge of naked DNA elastic con- 
stants, and to investigate quantitatively how they vary 
with the relaxed fiber structure, itself controlled by the 
"microscopic" structural parameters, as linker length 
or equivalently its relaxed twist r°, and entry/exit angle 
$ of linker DNA on the nucleosome . The results sub- 
stantiate an effective continuous description of the 30- nm 
fiber as an extensible worm-like rope (EWLR) extending 
the classical worm-like chain model to an extensible 
and twistable chain |15| |16| In Section V, numeri- 
cal results are presented and their physical meaning and 
implications are discussed. We then enlight our analyti- 
cal results by comparing them with experimental results 
obtained by Cui and Bustamante |^ by pulling a single 
chromatin fiber. Section VI presents the biological rel- 
evance of our study. In particular, we exploit our com- 
plete and quantitative analysis to discuss how the elastic 
properties of the chromatin fiber might at the same time 
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favour DNA compaction into the chromosomes and al- 
low local decondensation of chromatin involved in gene 
expression. 

II. ASSEMBLY OF A MODEL CHROMATIN 
FIBER 

Chromatin is composed of a double-stranded DNA 
molecule wrapped from place to place around histone 
cores We here focus on the fixed-nucleosome case, 
corresponding to deacetylated histone tails This al- 
lows to model separately the nucleosomes and the linkers, 
i.e. the naked-DNA segments connecting the successive 
nucleosomes. Our model amounts to build step-by-step, 
i.e. nucleosome after nucleosome, a chromatin fiber. It 
incorporates aknowledged data on DNA structure, its 
bend and twist persistence lengths, and on the nucleo- 
some structure. Each step of the assembly is described 
analytically, but we implemented this model within a 
Maple program in order to handle an arbitrary number 
of nucleosomes and to perform a quantitative analysis of 
the SH geometric properties. 

A. Linker modeling 

Linker DNA is in the usual B-DNA form: it is a right- 
handed double-helix, of radius awl nm and pitch equal 
to /° — 3.4 nm , which corresponds to — 10.6 base 
pairs (bp) per turn. We shall suppose in the following 
that all the linkers have the same number ribp of base 
pairs, i.e. the same length I = n^pl^ /n^^. It corresponds 
in vivo to phased nucleosomes observed for instance in 
satellite chromatine . It corresponds in vitro to recon- 
stituted fibers, involving repeated sequences each having 
a strong positioning effect on nucleosomes In 
fact, we shall need this hypothesis to be satisfied only lo- 
cally, over a few linkers, so that our model also accounts 
for irregular native chromatin. 

Since we are looking for generic properties of the chro- 
matin fiber, originating from its assembly, sequence ef- 
fects are ignored (besides, they might be treated in a 
second stage, within the same modeling, as local strains 
superimposed to the relaxed homogeneous structure that 
we here consider). Without sequence effects, the linker 
is straight in absence of applied constraints; experimen- 
tal evidence for straight linkers in relaxed fibers supports 
our generic modeling pl[ |. 

Linker DNA will be considered as an homogeneous 
cylindrical rope of radius a, that may twist and bend, but 
not stretch: due to the large value of the stretch modulus 
IDNA ~ 1200 pN 1^ (22|, stretch energy and strain can 
be ignored in linker DNA, at least in the low-stress situ- 
ations that will be considered here (forces below 20 pN). 
Linker DNA is thus seen as a non extensible semi-flexible 
polymer and described within the continuous worm-like 



rope model (worm-like chain model supplemented with 
twist energy [^); the elastic energy densities (energies 
per unit length) thus write: 

etw^st =^£^{lj-ljY C-75nm (twist) 

ebend = pIna a = 53 nm (bend) 

(1) 

where T is the temperature, ks the Boltzmann constant, 
Pdna the local curvature of the constrained linker and 
bj the twist rate. The relaxed twist rate lo^ = 27r/^° is 
supposed to be homogeneous, which again amounts to 
ignore sequence effects. The constants A and C are re- 
spectively the bend and twist persistence lengths of DNA, 
here given for lOmM NaCl (or any other monovalent 
salt) [0 [Q. We expect C to be almost indepen- 
dent of the ionic strength, since twist constraints are not 
strongly affected by electrostatics. By contrast, the elec- 
trostatic contribution to A is important ||2j] The 
non-electrostatic contribution provides a lower bound 
A > 23 nm, which is in fact an overestimate; the ex- 
perimental lower bound is rather A> AQ nm [ p6[ . 

We choose a direction along the DNA path correspond- 
ing to increasing indices j. We denote s G [0,/] the ar- 
clength along the dsDNA axis. Geometrically, the linker 
i is described as represented on Figure 1 by the local axis 
Uj{s) of the double helix and a vector tj(s), locally tan- 
gent to the minor groove and rotating around Uj{s) with 
an angular rate uj. tj{s) accounts for the torsional state 
of the linker: the angle of the rotation transforming t(0) 
into t(l) is precisely the twist t = Zcj of the linker. The 
vectors Uj{s) and tj(s) are normalized to 1 and make a 
constant angle: Uj{s).tj{s) = si'a{aDNA) independent of 
j and s. The slope of the strands in the double helix 
writes tan(Q;DAr^) — Z°/27ra, hence aoNA ~ 28°. The 
vector Uj{s) is independent of s in the relaxed fiber (a 
variation with s would describe the linker bending in con- 
strained state), whereas tj(s) rotates around itj with the 
relaxed twist rate 



B. Nucleosome modeling and assembly of the 
chromatin fiber 

Nucleosome structure is now well-known, thanks to 
high resolution cristallography data ||l^. Histone tails 
lock DNA on the histone core, so that nucleosomes are 
fixed on the DNA and the net effect of the nucleosome 
on the DNA trail can be described as a rigid kink con- 
necting the linkers j and j -|- 1. In consequence, twisting 
the linker j of a given angle reflects in a rotation of the 
same angle of the linker j -I- 1 around the linker j. The 
rotational positioning of nucleosome j + 1 with respect 
to nucleosome j is thus entirely prescribed by the twist 
angle t = Ico. In the relaxed state, it is equivalently 
prescribed by the linker length Z (cj = o;*^, t° = lui^). 
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Technicalities about the assembly and its geometric 
description are presented on Figure 2. The linker j en- 
ters the nucleosome j at an entry point Ej well-localized 
on the histone core. The linker DNA is grafted to the 
nucleosome by its minor groove ]l3t ; this implies that 
the plane spanned by Uj and tj{l) prescribes the position 
and the orientation of the nucleosome, given its radius 
fnuci = 5.5 nm, its height H = 5.7 nm and the angle 
a = 4.47° between the linker and the nucleosome axis 
Nj {uj.N = sin a). As shown in Figure 2 (right), uj, 
tj{l) and Nj belong to the same plane, which is nothing 
but the tangent plane to the nucleosome at Ej . Similarly, 
the plane tangent to the nucleosome j at the exit point 
Sj where linker j -I- 1 leaves the nucleosome is spanned 
by Uj+i and ij+i(0), which thoroughly determines this 
outgoing linker. The relative positioning of the linkers is 
conveniently described by the angle $ ($ > by conven- 
tion) between their projections on a plane perpendicular 
to the nucleosome axis Nj. The structural effect of nu- 
cleosome j amounts to a translation in space from the 
entry point Ej to the exit point Sj and a rotation of an 
angle 47r — $ around the nucleosome axis Nj , transform- 
ing {uj,tj{l)) into {uj+i,tj+i(0)). As long as we consider 
an homogeneous fiber, $ is independent of j. The line 
passing through the nucleosome center Gj and directed 
along the bissector Dj of this angle $ is the so-called 
dyad axis of the nucleosome j. Experimental evidence 
that the entry-exit angle $ varies with salt concentra- 
tion (it contracts when salt concentration increases) ||l9| ] 
pl[ [ p7[ motivates us to study various values of $. 
We here ignore interactions between nucleosomes. This 
is valid in low-salt conditions, in which the chromatin is 
expected to be extended, with internucleosomal distances 
larger than the interaction range. In high-salt conditions, 
our model aims at describing the "geometric" contribu- 
tion to chromatin elastic behavior, i.e. the contribution 
originating from the DNA elastic properties and relayed 
by the chromatin fiber architecture. Comparison with 
observed elasticity would give access to the neglected in- 
teractions. 



C. Linker histones and chromatosome modeling 

In vivo and in vitro experiments show that the con- 
densed conformations of the chromatin fiber are obtained 
for a more complex assembly of DNA and histones: an 
additional histone HI (or a close chemical variant H5, 
for instance in chicken erythrocytes) is bound to linker 
DNA near its entry/exit site on the nucleosome |2l| [ p9| . 
Presumably, the role of this "linker histone" is both to 
secure the wrapping of DNA around the nucleosome and 
to stabilize the DNA helical coiling forming the 30-nm 
fiber |3^]. Digestion experiments evidenced that the 
"core particle" now involves 166-168 bp of DNA (among 
which the 146 bp wrapped around the histone octamer); 
this new, larger, entity is called a chromatosome. The 



position of linker histone with respect to linker DNA and 
nucleosome, as well as the wrapping of DNA around it, 
are still debated |^ It is only recognized that 166- 
168 bp are protected and that the angle (still denoted 
$) between ingoing and outgoing linker DNA is lowered, 
varying with ionic strength. Our model easily accomo- 
dates the presence of linker histones along the fiber since 
only the resulting kink of DNA path is to be described. 
We may roughly account for the presence of linker his- 
tone by modifying $ = 130° (without linker histone) into 
$ = 90° at low salt down to $ = 45° in high salt |2^] 
plf ; the nucleosome is then replaced, within the same 
two-angle modeling, by an effective cylindrical chromato- 
some, with possibly different values for Vnuci, H and a. 
For instance r„uci = and a lower value of H would re- 
produce the crossing of ingoing and outgoing linkers in 
the neighborhood of the linker histone. Also, the effective 
length of the linkers is to be reduced by 20-22 bp. 

We might thus study within the same geometrical mod- 
eling all the different instances encountered for the nucle- 
osome: basic nucleosome or nucleosome dressed with HI, 
with or without chemical modification of HI tails. The 
specific biological details (HI positioning, possible acety- 
lation of histone tails |^^) would be taken into account 
first in whose value is experimentally accessible, then 
in a second step by a precise fit of the parameters r„„c;, 
H and a of the effective cylindrical "core particle" , so 
as to reproduce accurately the kink induced in the linker 
DNA broken line. Nevertheless, a systematic study led 
us to think that $ is a far more important control pa- 
rameter than Vnuch H and a, which will be henceforth 
considered as fixed characteristics. 

To summarize, the parameters of the model, to be var- 
ied in our study, are the two angles (or equivalently 
I = /lo^ or nbp) and $. The angle determines the rel- 
ative positioning of two successive nucleosomes, whereas 
the angle $ determines the relative positioning of two 
successive linkers. 



III. THE CHROMATIN SUPERHELIX (SH) 

The regularity of the microscopic assembly enforces an 
helical organization: we shall speak of the chromatin su- 
perhelix (SH). We now describe the main geometric prop- 
erties of this helical coiling resulting from the regular as- 
sembly of linkers and nucleosomes. As shown in Figure 
3, our qualitative structural predictions recover standard 
structures proposed yet long ago |Q Q. In particu- 
lar, we note that the various structures obtained when 
varying r° (or equivalently ni,p) and $ cover both ribbon 
structures (less than 3 nucleosomes per turn) and cross- 
linked structures (more than 3 nucleosomes per turn). 
We avoid to speak of "solenoid" structure as this term 
refers to a model with bent linkers , whereas our mod- 
eling considers that linkers are straight in the relaxed 
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fiber. 



A. Symmetry properties of the SH 



described by the angle /? G [0, 27r[ between its axis Nj 
and the SH axis A, in the frame {A,D,A A D). This 
angle satisfies cos /3 = Nj .A. 



We first underline the symmetry properties of the SH, 
of much importance since they will reflect in the elas- 
tic properties of the fiber in the linear regime. The SH 
exhibits a discrete helical invariance: all the linkers (re- 
spectively all the nucleosomcs) are equivalent, up to a 
rotation of angle 6 around the axis A of the SH and a 
translation along this axis, bringing the linker j onto the 
linker j + \ (respectively the nucleosome j onto the nu- 
cleosome j + 1). The nucleosome centers all lay on a 
cylinder of axis A and radius R, whereas the entry and 
exit points all lay on a cylinder of radius r (see Figure 
4). 

As the two strands of dsDNA are oriented in opposite 
directions (see Figure 1), there is no privileged orienta- 
tion along linker DNA and only the nucleosome geome- 
try can break the symmetry upon reversal of the fiber. 
Hence the fiber is not oriented provided the entry and exit 
points in a nucleosome have the same characteristics: it 
is the case in absence of linker histonc, or when linker 
histone is positioned symetrically with respect to Ej and 
Sj. The fiber properties are then unchanged under the 
transformation: 
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(2) 



Here j = . . . N + 1 where N is the number of nucle- 
osomes in the fiber beginning in 5*0 with linker 1 and 
ending with linker A''-!- 1, at En+i It follows that the nu- 
cleosome axes are all tangent to the cylinder containing 
the nucleosome centers, and that the lino Dj relating the 
nucleosome center Gj to its orthogonal projection onto 
the SH axis is nothing but the dyad axis of the nucleo- 
some j (see Figure 5, left). This dyad axis is invariant 
under the above reversal transformation. When r° is 
varied, the linker direction Uj sweeps a cone of axis Dj 
(or -Dj-i, due to the symmetry upon reversal of the SH) 
and angle ^, depending on $ (see Figure 5, right). This 
property simply reflects the fact that the geometry of the 
fiber assembly is thoroughly determined once the nucle- 
osome is positioned with respect to the SH axis, due to 
the symmetry properties and to the fact that the junc- 
tion between the nucleosome and the ingoing and out- 
going linkers is a rigid kink. Indeed, by symmetry, the 
nucleosome dyad axis Dj is orthogonal to the SH axis A 
and points towards it; the remaining "degree of freedom" 
is the orientation of the nucleosome axis Nj in the plane 
orthogonal to Dj ; when I varies, Nj rotates together with 
Uj and Uj + 1 around Dj. This orientation is thoroughly 



B. SH geometric characteristics 

The relevant quantities to be computed (lengths D and 
r, angles 0, (3, z, rj, and ^) are shown on Figures 4, 5 and 
6. All these quantities are independent of j, due to the 
(discrete) rotational invariance of the SH. When linker 
lengths are all increased by the DNA pitch the twist 
angle of each linker increases by 2n, hence the relative ori- 
entations of successive linkers, or successive nucleosomes, 
are unchanged. It follows that the SH shape is preserved: 
the angles (as z, r] or 6) are unchanged whereas the length 
r is scaled by a factor (1 -|- l^/l) (when comparing I and 
l + l'^); the lengths Rsh and D are also scaled, but the in- 
volved scaling factor is smaller, since the contribution of 
the nucleosomes to the actual size of the SH is the same 
for I and I + l'^; the precise value of this factor moreover 
depends on the structure of the SH (nucleosome orienta- 
tion, for instance). 

We here underline two points. First, the orientation of 
the SH axis is chosen so that the distance D between two 
successive nucleosomes along the SH axis A is positive: 
D > (i.e. D = A.[GjGj+i] > by definition of A). 
Also, since the structure is discrete, the angle 6 of the 
^Ttitation TZg around the SH axis A, transforming the pro- 

Kction of a nucleosome onto the plane orthogonal to A 
following one is not univoquely defined: 6*0 > 
and ^0 — 27r < are both possible (think to the array of 
grains in a corn ear or to the cells on a pineapple peel). 
By convention, we choose the value of smallest modulus, 
hence let vary in ] — tt, tt]; it corresponds to the (pro- 
jected) angle swept by the geodesic (the array of grains 
of largest slope) relating Gj to Gj+i onto the cylinder of 
axis A and radius R. 

A= dist(Gi,G2) and the angle * = G^G^s defined 
by three The pitch P of the SH is given by P = 27rZ)/|6»|. 
Note that the choice of the value of smallest modulus for 
9 is essential to get the actual value of the pitch through 
this formula. The degree of compaction can be measured 
as the number of nucleosomcs per 10 nm of fiber, equal 
to 10/ D when D is given in nm. 



C. Excluded- volume effects 

Excluded volume (of nucleosomes but also linkers) has 
to be taken into account to discard unrealistic struc- 
tures. In consequence, the above results should be sup- 
plemented by a geometric condition assessing whether 
they actually correspond to a possible structure, i.e. 
compatible with excluded-volume constraints. We here 
rather state an upper bound, involving only the geomet- 
ric parameters of the fiber. Recall that /3 is the angle 
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between the nucleosome axes Nj and the SH axis (inde- 
pendent of j by symmetry). A sufficient condition en- 
suring that no steric hindrance is encountered in the SH 
structure is that P > Pc where: 

= i^cos/3 + 2i?|sin/3| (3) 

where R is the distance of the nucleosome centers to the 
SH axis (computed in the Maple implementation) and H 
the height of a nucleosome. For /3 = 0, this upper bound 
gives the exact threshold Pc = H . All linker lengths 
such that P{1, $) < Pc are forbidden. Note that a lo- 
cal breaking of this criterion is allowed: it is still possi- 
ble that a few linkers lengths (less than 2-1:/ 9) take this 
value. A more accurate check can be obtained by draw- 
ing for each configuration the developped representation 
of the SH, i.e. the mapping on a plane of the fiber track 
on the cylinder of radius R, as shown on Figure 7. Its 
construction is straightforward, knowing the geometric 
characteristics D, 9 and (3. Excluded-volume constraints 
are satisfied if the tracks of any two nucleosomes do not 
intersect. This method requires to check each configura- 
tion, hence comes after the criterion on P ~ Pc, when it 
fails. 



D. Numerical results 

The major interest of our numerical implementation is 
to go further than a qualitative description and to com- 
pute explicitly any geometric characteristic of the relaxed 
fiber, in order to study quantitatively its variation with 
I (or r'^) and $. The main results are presented in Fig- 
ures 8, 9 and 10 for a linker length between 30 and 50 
bp (usual case). The quantitative geometric character- 
ization of the SH is essential to compare the different 
structures obtained for different values of (Z,$). It is 
also interesting from an experimental viewpoint, since 
it would allow to extract informations on the chromatin 
assembly from the experimentally accessible values of D 
and 9. 

We first evidence that a regular helical packing of nu- 
cleosomes is compatible with straight linkers and does not 
require nucleosome-nucleosome interactions (but the spe- 
cial configurations exhibiting columnar arrays of nucleo- 
somes are likely to be stabilized by such internucleosomal 
interactions). Mainly, our systematic analysis shows that 
it always lead to a 30-nm fiber. Although the detailed 
structure of the fiber is strongly sensitive to linker length 
I and entry/exit angle $, the value of about 30 nm for 
the fiber diameter 2Rsh is a robust feature, hence of low 
structural significance; for instance, it does not discrimi- 
nate close cross-linked structures and more extended rib- 
bons; 

We enlight some remarkable geometries (see also Sec- 
tion V.E): 

• a ribbon - or zigzag - structure is obtained for 9 ~ n, 
when the number of nucleosomes per SH turn reaches 



the minimal value 2; then all the nucleosome dyad axes 
are aligned, P is close to tt/2 but slightly different: /3 = 
tt/2 ~ a (see Figure 14). 

• a particular packing is obtained when THip/n^^ equals 
an integer; the nucleosome axes are then all aligned, 
hence all parallel to the SH axis: /3 = 0. A fine tuning of 
riiip or $ brings then into columns as seen on Figure 3. 

• the vanishing of the pitch P is accompanied by a rever- 
sal of the SH axis, from which follows that cos z exhibits 
a jump from a value c to the value — c and 9 jumps to 
the value 2ti ~ 9. We shall see that the vanishing of the 
pitch P is an important (although virtual) event, having 
striking implications on the elastic properties of the fiber. 

• a change of chirality is observed for 77 = tt: the SH is 
right-handed for r] G [0, tt] and left-handed for 77 G [vr, 27r]. 
For r/ = TT, linkers cross the SH axis. 

One has to carefully distinguish between: 

- the handedness of the track passing through the nucle- 
osome centers. As we have seen above in the definition of 
9, the handedness of this discrete structure is ill-defined. 
A same architecture can be seen as right-handed (direct 
rotation TZg) or left-handed (indirect rotation 'R-g-2Tr)- 
The choice of the value of smallest modulus, adopted 
here, corresponds to the shortest path. This "pseudo- 
chirality" is then given by the sign of 9 and changes 
when 9 crosses the value ±7r or when D (or cquivalcntly 
the pitch P) vanishes. We underline that this structural 
feature does not define a relevant chirality, due to the ar- 
bitrariness of its definition, and should not be confused 
with the chirality of the SH, defined as follows; 

- the chirality of the linker DNA trail (a broken line) ; it is 
determined by the position of 77 with respect to tt, i.e. by 
the sign of cos(?7/2). As fiber elasticity originates from 
linker elasticity, this chirality is the only relevant one 
for our mechanical study and will be henceforth adopted 
as the definition of the fiber chirality. We expect this 
chirality to determine the coupling between the twist and 
stretch elastic degrees of freedom i.e. the sign of g to 
coincide with the chirality. 

In the same spirit, we warn about the difference be- 
tween the fiber and an helical DNA coiling passing 
through the nucleosome centers, in particular in what 
concerns their mechanical properties. 

One of the conclusions of this thorough structural 
study is the fact that the connection between the mi- 
croscopic parameters and the SH geometric characteris- 
tics is too complex and multivariate to get small-scale 
informations from structural observations of the fiber. 
We underline that the study of chromatin structure is 
not in itself sufficient to unravel its biological function- 
ing, all the more as structural observation of the fiber 
is difficult and some results questionable |^^. Moreover, 
structural changes might not be the only way to pass lo- 
cal DNA modifications to higher scales; the functioning 
might rather be controlled by elastic properties. 

For these two reasons, we turn to the analysis of the 
mechanical properties of the fiber, which are now at hand 
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since we know explicitly all the geometric characteristics 
of the fiber assembly. 



IV. ANALYTIC CALCULATION OF THE 
ELASTIC CONSTANTS OF THE CHROMATIN 
FIBER 

A. Modeling the chromatin fiber as an extensible 
worm-like rope (EWLR) 

Our aim is to study the elastic response of the fiber to 
external stresses at scales larger than its pitch. We thus 
consider that we apply to the fiber a force F, along the 
fiber axis A; if F were not directed along the SH axis, it 
would induce a transitory motion, compelling the axis to 
align itself with the force direction, at least far from the 
ends. We here only suppose that this constrained equi- 
librium is yet reached. We may also apply a torsional 
torque Mt (directed along A) and a flexural torque Mb, 
i.e. a torque component orthogonal to A. Computation 
of elastic coefficients describing the linear response of the 
fiber to the applied force and torque will be performed 
analytically within a continuous description of the fiber, 
i.e. an effective large scale description in which the dis- 
crete nature of the assembly is smoothed out. We denote 
S the arclength along the axis of the relaxed SH, u{S) 
its local relative extension, ^1{S) its local twist rate and 
g{S) its local curvature, with respect to the straight and 
untwisted relaxed state of the fiber for which u, Q and 
p thus identically vanish. Due to the axial symmetry of 
the relaxed SH at scales larger than its pitch, the bend- 
ing energy density does not depend on the direction of 
the bending so that it involves only the total curvature 
p and the modulus A/f, of the fiexural torque (the spe- 
cial instance of ribbon-like configurations that break this 
axial symmetry is discussed below in Section V-E). The 
fiber has thus only three degrees of freedom u, Q and p, 
local along the fiber. These strains u, fl and p are the 
canonical variables of the density of elastic free energy, 
which reflects in the following differential form: 



desH = F du + Mt dVl + Mb dp 



(4) 



We restrict to the linear response regime, which expresses 
in the following relation: 
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where it can be shown that the stress-strain tensor T 
is necessarily symmetric (a special instance of Onsager 
relations). This linear response ansatz relates the strains 
u, and p of the fiber and the stresses F, Mt and Mb 
experienced by the fiber. A is the bend persistence length 



of the fiber, C its twist persistence length, 7 its stretch 
modulus (dimension of a force) and g the twist-stretch 
coupling constant (no dimension). 

Plugging the linear response ansatz (H) into (^ leads 
to the SH density of elastic free energy: 



keTAg^iS) keTCn^ lu^{S) 



kBTgn{S)u{S) 
(6) 



Such a continuous description of the fiber can be termed 
"extensible worm-like rope" model (EWLR) It ex- 
tends the WLC model introduced in 1949 by Kratky and 
Porod and currently used to describe stiff polymers JM] 
by accounting for twist (as above for linker DNA [p,2||) 
but also stretch degrees of freedom. The fact that chro- 
matin is chiral demands a linear coupling between twist- 
ing and stretching. Due to the axial symmetry of the SH, 
{F,Mt,Mb) and {F,Mt,-Mb) (and also {F,Mt,nMb) 
where TZ is any rotation around the axis A leaving the 
SH unchanged) should induce the same energy change 
in the SH. This symmetry argument shows that there is 
no other coupling term at the linear order considered 
here This model has been fully investigated by 

many groups in the context of DNA; in particular, force- 
extension curves have been obtained p^]; a nonlinear 
term V{u) might be added in (||) to go beyond the linear 
response regime. These results can be straightforwardly 
transposed to chromatin fiber, so that the description of 
the harmonic elastic behavior of the chromatin fiber re- 
duces to the computation of the four elastic coefficients 
A, C, J and g involved in the EWLR model. Their de- 
termination from the computation of the elastic energy 
stored in the constrained linkers will prove that the fiber 
actually fits in an EWLR model when considered at large 
enough scale. 



B. A general analytic method for computing elastic 
coefficients 

The first aim of our study is to express the elastic co- 
efficients A, C, 'y and g as a function of the elastic coeffi- 
cients of linker DNA, given the relaxed geometry of the 
fiber. A key point of our approach is to relate the stresses 
exerted respectively on the linker, considered as a WLR, 
and on the fiber, considered as an EWLR instead of re- 
lating the strains {uj — lo^ , pjjj^a) of the linkers and the 
strains (u, f2, p) arising at the fiber level (this approach 
appears to be technically cumbersome and analytically 
intractable, see Section VII. B). Although the stresses F, 
Mt and Mb are not the canonical variables of the free 
energy density esH^ we shall express it as a function of 
these variables F, Mt and Mb- Plugging into (^) leads 
to the following expression of the free energy density: 



[ksTCF'^ 



77\f2 - 2kBTgFMt\ 



M^b 



2{kBTCj - klT^g^) 



2kBTA 



(7) 
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The principle of the computation is the following: know- 
ing the relaxed geometry of the fiber, i.e. the quanti- 
ties r, ri, z and D computed in Section III, it is possible 
to analytically determine the local stresses experienced 
by linker DNA when the fiber is constrained at its ends 
(and only at its ends) and to deduce the elastic energy 
of a linker as a quadratic function of {F, Mt, Mi,). The 
identity of the elastic energy of the fiber and the elastic 
energy stored in the linkers will allow a termwise identi- 
fication with (0), which leads the values oi A, C, 7 and 
g as analytical functions of the geometric parameters of 
the fiber, and through our numerical implementation, as 
a function of / (or r°) and 

A much important intermediate result is to establish 
the relation between the local stresses (experienced by 
the linkers) and the global stresses (applied on the chro- 
matin fiber). The current point on linker j, with ar- 
clength s, is denoted Pj{s). Recall that linker j leaves 
the nucleosome j — 1 at Sj-i and enters the nucleosome 
j at Ej (see Figure 6). We denote fj{s) the force and 
TOj(s) the torque exerted at the point Pj{s) of linker j 
by the upstream part of the fiber. The main arguments 
involved in the derivation are the following: 

- when only pulling the fiber, the rotational invariance 
of the relaxed fiber should be preserved when end effects 
are ignored. We may (and in fact should) focus on the 
universal, rotationally symmetric behaviour of the fiber, 
observed far from its ends. Indeed, the quadratic, rota- 
tionally invariant EWLR energy is to be fitted only to 
the quadratic, rotationally invariant energetic contribu- 
tion coming from the linkers, otherwise the identification 
would not make sense. 

- we restrict to the linear regime. We may consider 
separatedly the different stresses applied to the fiber and 
simply sum up their effects to recover the general solu- 
tion. The contribution of the applied torque M to the 
local torque 'rhj{s) is merely M. 

The conclusion follows using standard equilibrium 
equations of spring mechanics. When external forces and 
torques are applied at the ends only, linear response hy- 
pothesis implies that the solution finally writes: 



m = F 

m,{s) = M-[0,{s)P,{s)]hF 



the fiber architecture in the expression of the local torque. 
As we restrict to the description of harmonic elasticity, 
the coefficients of F and M in the elastic energies of the 
linker will be computed within the relaxed SH. We un- 
derline that we do not need to compute the constrained 
shape of the linkers to describe the linear response of the 
fiber to applied force and torque. We carry on the compu- 
tation by exploiting the fact that the fiber is not oriented. 
It follows that: \\0{E^j)Ej\ \ = 110(5^)5^11 = r and: 

\\[0,{s)P,{s)]AA\f = \\0,{s)P,{s)\f^r'[cos'{rj/2)co8'z + sm^{Tj, 

(9) 

where 77, r and z are the values associated to the relaxed 
SH. 



C. Small-scale grounds for the EWLR modeling of 
the chromatin fiber 

Linear response ansatz applies both at the linker level 
(WLR model, with parameters A and C) and at the SH 
level (EWLR model with parameters A, C, 7 and g). 
Knowing the stresses thus gives the strains, respectively 
(pdna, w — uj°) and (m, ft, p), if required, and the elastic 
energies. The computation of the SH elastic constants 
rests on the obvious but essential fact that the elastic en- 
ergy stored in the SH, expressed by (@) at the SH level, 
is nothing but the sum of the elastic energies stored in 
its linkers, since no interactions are involved in our mod- 
eling. 

A key point is the different decompositions of the 
torques into torsional and flexural components at the 
fiber level and at the DNA level: 



M{S) = MtA + Mb{S) (SH) 



m(s) — mtu + : 
(10) 



(8) 



The relation between the two decompositions involves the 
fiber geometry. We underline once more that in the linear 
response regime, only the relaxed geometry is involved. 
The first decomposition gives the components involved in 
the EWLR energy density whereas the second is required 
to compute the elastic energy of a linker. Indeed, the 
elastic energy densities of a linker can be expressed as a 
fimction of local stresses, according to: 



where Oj (s) denotes the orthogonal projection of Pj (s) 
onto the SH axis A. It gives the relation between the 
global stresses F and M exerted at the fiber ends and 
the local stresses experienced at the linker level, at each 
point of the DNA path. The rotational symmetry of the 
fiber ensures that F is directed along the SH axis A at 
equilibrium. 

The term [Oj{s)Pj{s)] A F refiects the involvement of 



mt = kBTC{uj — bj^) 
nib — ksTA pDNA 



2kBTC 



and efce 



(11) 



According to a general result ||], besides easy to check 
directly here, mt is constant along the linker (i.e. in- 
dependent of s) whereas rfib may vary with s both in 
direction and modulus. 
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The EWLR model is a continuous model. It makes 
sense to describe the elastic behavior of the discrete 
chromatin structure by means of an EWLR model only 
at large enough scale, so that the discrete effects are 
smoothed out. At lower scale, the specific orientation 
of each linker influences its energy, more precisely the 
contribution coming from the global stress Mf,, which 
breaks the rotational invariance; the persistence length 
of bending A of the SH is only defined after averaging 
over, say, one turn of SH. At the linker scale, twist-bend 
and stretch-bend coupling terms are present; they vanish 
on the average. The averaging keeps only the large scale, 
rotationally invariant, "EWLR-like" contribution. The 
remaining terms describe local contributions to the bend- 
ing energy of the SH, cancelling each other, hence with 
no resulting effect at the fiber scale. The average over 
a number of linkers sufficient to recover the rotational 
invariance of the fiber (using J2f=i ^] ~ ('^^^ ^)^) 
conveniently replaced (except for the ribbon-like struc- 
ture for which the rotational invariance breaks down, see 
Section V-E) in the computation by an average over the 
directions of Mb (denoted <>): it has the same effect of 
extracting only the resulting contribution at large scale. 
After integration of the densities tbend and Ctwist along 
the linker, the elastic energies Eb(,nd and Etwist stored in 
a linker write: 



A 



C = 



AD /I 



CD ( 



tan^(77/2) 
3 



cos^ z + 4; sin^ z 









3v 



kBT 



A cos z-\-C sin 
r'^ cos^(r//2) 



-A) sin z cos z 
r cos(r^/2) 



H- 



tan2(^/2) 
3 



(16) 



We underline that our approach yields the elastic con- 
stants of any relaxed geometry, hence allows an analysis 
of sensitivity with respect to / (i.e. r") and <&. 



V. RESULTS AND PHYSICAL DISCUSSION 



A. Numerical results for the elastic constants 



1 



2AkBT 

GAkeT 
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2CkBT 
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< ml{s) > ds 

Hi) 



(12) 



3 cos zcos' 
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< TO, > ds — 
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mt > 



2CkBT 



2CkBT 



(^Fr cos (J^^ sin z -I- Af, 



t cos z 



Numerical implementation of the above analytical for- 
.^.^j .n^las cffiAfigjperfornied for fiajy\valuefeMf the micros^om 
2AkBT ^'f?ai%,ihe'^jf^fia^,^&Sffi]lt^igij;idii^Jg^6A l'rf;93iaM;s 
to a Maple program. We here present the results ob- 
tained in the two situations $ = 90° and <i> — 50", taking (14) 
as a variable the linker length, 
sm' z 

Starting from the microscopic structural parameters (15) 
and <i>, passing through the determination of the SH 
geometric characteristics D, r, z and rj, we obtain the 



ACkeT 



D. Analytic expression of the elastic constants as a 
function of geometric parameters of the fiber (SH) 

The quadratic expression of Ebend + Etwist as a func- 
tion of F, Mt and Mb is straightforwardly identified with 
(^ integrated over the length D. This justifies a poste- 
riori to map the fiber and its elastic behaviour onto the 
continuous EWLR model. In other words, it provides 
a microscopic validation of the EWLR modeling of the 
chromatin fiber. In particular, we check the expected 
vanishing of the twist-bend and stretch-bend coupling. 
Some algebra finally yields the elastic constants of the 
SH fiber: 



explicit values of the elastic constants, as functions of r*^ 
(or equivalently Ubp) at fixed (arbitrary) $. Indeed, ex- 
pressions (^6|) are valid for any relaxed structure hence 
allow to analyze the complete range of variations of the 
elastic constants as r° and $ vary, as shown in Figure 1 1 . 
We here present two typical cases, currently proposed in 
chromatin structural studies: $ = 90° then $ = 50°. A 
decrease of $ is presumably induced by an increase of salt 
concentration, lowering the mutual repulsion of ingoing 
and outgoing linkers and strengthening the structural ef- 
fect of linker histones. Experimental observations of a 
fiber dressed with linker histones showed that a value of 
<& around 90° corresponds to a salt concentration of less 
than 5mM NaCl, whereas <& decreases below 50° above 
15 niM NaCl ^ Q. 

The curves of Figure 11 evidence a strong sensitivity 
of elastic coefficients with respect to the fiber structure, 
as controlled by ribp. A striking result of our study is the 
sharp decrease of all elastic coefficients together with the 
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pitch P, around a critical value ntp^d^) for which P van- 
ishes. The critical value Ubp^d^) depends on the value 
of $ and more generally on the precise modeling of the 
chomatosome. As P vanishes, steric hindrance between 
nucleosomes precludes to build a regular structure with 
n-bp — nbp^d^), so that this feature might be seen as irrel- 
evant. Nevertheless, it is possible to have Ubp = Ubp^d^) 
locally, i.e. over less than one turn of SH; this flexibility, 
although local, might yet have dramatic consequences: it 
is possible to create a noticeable kink in the chromatine 
fiber at this point by applying only very weak stresses. 
Such "critical" turns appear as defects where the regular 
compact fiber is easy to "open". Moreover, the curves 
near Ubp^c but above the excluded-volume threshold Pc 
are still influenced by their vanishing in Ubp^c- The steep 
variation of the elastic constants in the neighborhood of 
Ubp^c is an actual property, corresponding to a highly sen- 
sitive fiber: the geometric parameter D and the elastic 
constants increase by a noticeable factor within one or 
two base pairs. Another important result, whose implica- 
tions are discussed below (Section V-D), is the observed 
change of chirality of the SH, defined in Section III.D 
as the sign of cos(77/2) and exactly correlated with the 
changes in the sign of the twist-stretch coupling g. 

Finally, the relative energetic contribution of linker 
bending and linker twisting are shown on Figure 12 in 
pure cases where only one of the stresses F, Mt or Mb 
is applied (stretch, twist or bend at the fiber level). We 
checked that the stretching energy and the twist-stretch 
coupling energy of the linker are negligible (almost a few 
percent of the total clastic energy stored in a linker) due 
to the high value of ^dna (around 1200 pN) and goNA 
(between 20 and 30) ||l^ [g2|. The figure evidences an 
interesting feature coming from the spatial organization 
of DNA into a chromatin fiber. According to the fiber 
structure (i.e. linker length / and angle $) and to the na- 
ture of external stresses (pulling force, torsional or flexu- 
ral torque) the linkers react either by bending, either by 
twisting, which might play a biological role by regulating 
the local DNA structure. 

Our computation of A is based on its "energetic" def- 
inition through the expression (||) of the free energy of 
the EWLR description of the fiber: the energetic cost 
required for bending the SH axis and creating a uni- 
form curvature p is kBTAp^/2 per unit length. It orig- 
inates from the cost of twist and bend distortions in- 
duced at lower scale on the linkers. This definition does 
not (and does not have to) consider possible steric hin- 
drance (hard-core interactions). Indeed excluded-volume 
constraints do not contribute to elastic energies hence do 
not modify A] they only forbid some deformations. To 
determine whether a given bending of the SH is possible, 
one has 

- first to check whether such a bending is "geometri- 
cally" allowed, i.e. satisfied excluded-volume constraints 
( a necessary condition is that pRex < 1, i.e. the radius 



of curvature 1/p should be greater that the excluded- 
volume radius Rsh WA i 

- secondly to check whether enough energy is available, 
using the value A determined above, notwithstanding the 
excluded-volume effects. 

Note that these criterions are independent (geometric 
and energetic respectively) and they can be checked in 
any order. 



B. Comparison with an helical coiling (ordinary 
spring) 

It is interesting to compare the SH elastic behavior 
to that of a simple helical "spring" (a toroidal coiling 
of DNA), i.e. the continuous structure characterized by 
the same angle z (here between its axis and the tangent 
to the elementary fiber) and the same radius r. The 
helical spring is thus implicitely controlled by the same 
parameters I and $ as the SH. The comparison enlights: 

- that the chromatin fiber is actually a spring from 
a mechanical point of view, although a special, tunable, 
one; 

- that the difference between ordinary springs and 
chromatin lies in the angle ry (77 = in ordinary springs 
whereas it is not small in chromatin), from which origi- 
nates the tunable character of the chromatin spring. 

The reasoning to compute the elastic constants of an 
helical spring is at each step analog to the above one 
(Section IV). We obtain the following expression for the 
average energies per unit length: 



[Fr sin z + Mt cos z) 

2CkBT 
cos^ z F^r^ 



2AkBT 



M^ sin^ z 
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Ml sin^ z 
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A AkBT 



where the average is performed either on the orientations 
of Mb^ either on a turn, i.e. over an arclength 27rr/sinz 
corresponding to the pitch 27rrcotz. Such an average 
is essential to fit in the definition of elastic constants 
(continuous, rotationally invariant EWLR model). The 
expression of the elastic constants follows: 



Ordinary spring 
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A cos z 
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(19) 
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In fact, these classical formulas of spring mechanics can 
be recovered directly from those obtained for the chro- 
matin SH by letting Z — > at fixed z in the SH formulas; 
accordingly, — *■ and Djl — > cosz. This link proves 
that chromatin is actually a special kind of spring. In 
the chromatin fiber, sin(77/2) is not small and can even 
reach 1. Hence all the features depending on r\ will ex- 
hibit a striking difference when comparing the SH and 
the helical coiling of DNA. This is the angle rj ^ 0, re- 
flecting the discrete character of the chromatin assembly, 
which is mainly responsible of the tunable elasticity and 
tunable chirality of the chromatin fiber, as can be seen by 
comparing the elastic constants of the helical spring and 
of the chromatin fiber. 

Let us also compare the distribution of DNA elas- 
tic energy between twist and bend degrees of freedom. 
The comparison is presented in Figure 12. When only 
a torque (either torsional, either flexural) is applied, the 
partitions in the SH case and in the case of an helical 
coiling are exactly identical. When only a pulling force 
is applied, the energy is mainly stored in the torsional 
energy in the case of a simple helical spring whereas it 
is the converse in the chromatin SH. In the chromatin 
fiber, the energy stored in the linker bending is in any 
cases larger (and sometimes far larger) than the energy 
stored in the twist. This remark reflects the following 
property of the SH: given mt and rub (modulus averaged 
along a linker), it is possible to determine a force F and 
a torsional torque Mt creating mt and rub on each linker 
only if TOf, > Cte.mt, where the constant depends on the 
SH geometric characteristics. In particular, it is impos- 
sible to have mb = unless mt = 0, i.e. F = Mt = 0. 

This difference between the elastic behavior of the SH 
and the corresponding simple helical coiling originates 
from the discrete structure of the SH, where the link- 
ers are spaced by rigid kinks (the nucleosomes). This 
strongly modifies the relation between the pulling force 
applied at the fiber ends and the local flexural (m;,) and 
torsional {mt) components determining the DNA elastic 
energy. The discrepancy is maximal for sin(?7/2) = 1, 
which corresponds to the SH conformation in which the 
linkers cross the SH axis and SH chirality changes, which 
obviously never occurs in an helical coiling of bounded 
radius(either left-handed, either right-handed, but never 
in between if r neither vanishes nor diverges). Although 
the SH " looks like" an helical coiling (see Figure 3), its 
response to a pulling force differs from that of an helical 
spring. Indeed, we again underline that its elastic prop- 
erties are determined by the characteristics of the broken 
line formed by the linkers, which might be far different 
from an helical coiling (especially when r] is close to tt, 
which corresponds to cross-linked structures). 

Elasticity theory of an homogeneous cylindrical rod H] 
predicts the relation AAksT = ^Rg^- The predicted 
value of A using this formula is roughly correct for the 
helical coiling, but it overestimates the actual value of 
the SH persistence length by a factor of 2 or more. This 
discrepancy enlights the importance of the complex sub- 



structure of the SH, leading to a fiber far different from 
an homogeneously filled cylinder of radius Rsh and even 
from an helical coiling of DNA. 

C. The origin of the sensitivity of the elastic 
constants with respect to the fiber structure 

Although the chromatin SH and an helical coiling of 
DNA exhibit quantitative differences, their elastic behav- 
iors share the same key feature (see Figures 11 and 13): 
a sharp decrease towards of the elastic constants, to- 
gether with the pitch. Considering a simple helical spring 
yields the - quite intuitive - explanation of this feature: 
when the slope of the elementary (DNA) fiber decreases 
to (DNA path almost orthogonal to the fiber axis), 
the number of turns per unit length along the fiber axis 
increases (see in particular the divergence of the com- 
paction ratio 10/-D(nm) on Figure 9, outside the figure 
frame). This accumulation of turns, each acting as an 
hinge, easily allows a huge deformation of the fiber, sim- 
ply by the addition of small changes occuring in each 
turn. This explanation, involving only an accordion-like 
behavior, ensures the robustness of the feature. The flex- 
ibility dramatically increases when the number of coils 
per unit length of fiber increases, whatever the precise 
geometry with which the elementary fiber is coiled into 
a large-scale one. 

The strong correlation between the pitch and the four 
elastic constants behaviors when ni,p varies simply re- 
fiects in the very expression of the elastic constants: they 
all write as (D/l) f{z, r), hence are slaved to the vanish- 
ing of D (or P). The factor D/l shows that the extensibil- 
ity and flexibility of the SH are controlled by its degree of 
compaction. The more compact the conformation (low 
value of P), the more flexible is the fiber, easily bent, 
twisted or stretched. We claim that the result is typical, 
i.e. does not depend on the details of the modeling, but 
only on the general feature of the chromatin fiber assem- 
bly, i.e. nucleosomes linked by segments of DNA. This 
form {D/l) f{z,r) also indicates that the elastic proper- 
ties of the fiber are mainly controlled by the angle z be- 
tween the linkers and the SH axis. Indeed, z roughly de- 
termines the decomposition of local stresses among twist 
and bend degrees of freedom of linker DNA (stretch is 
always negligible). 

D. A multi-strand spring with tunable chirality 

We first underline that the chromatin fiber is a spring 
of tunable chirality: the sign of the twist-stretch cou- 
pling g, associated to the chiral nature of the SH, changes 
twice within a "period" (an /-interval of width 10.6 bp, 
i.e. the DNA pitch), see Figure 11; correspondingly, a 
change from a right-handed SH to a left-handed SH oc- 
curs, as seen on Figure 1. Reverse behaviors are observed 
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in right-handed and left-handed SH. Indeed, change of 
chirahty exactly corresponds to change in the sign of the 
susceptibility dD/dr, i.e. the slope of the curve giving 
D as a function of Ubp (see Figure 8), up to a factor 
dnhp/dr = n^^/^'K. The noticeable consequence is that 
a change At will either decondense or compact the fiber 
according to the sign of dD/dr, i.e. to the chirality of 
the SH. 

Far more, the two- level structure of the chromatin fiber 
makes it quite similar to a multi-strand spring: an in- 
teresting behavior arises from the interplay between the 
right-handedness of the linker DNA and the tunable chi- 
rality of its coiling within the SH. Applying (only) a 
pulling force FA creates a torsional torque: 



•nit = Fr sin zcos{ri/2) 



(20) 



The sign of cos(r//2) thus determines the sign of (re- 
call that by definition sinz > 0). According to this sign, 
i.e. to the SH chirality, pulling the chromatin fiber will 
either twist, either unwind the linker DNA double he- 
lix. This behavior is well-known in the context of multi- 
strand springs . Its biological interest might be to reg- 
ulate DNA dcnaturation required for transcription and 
replication. We get a similar result when applying a tor- 
sional torque MfA; the torsional torque at the linker level 
writes rrit — Mjcosz. A change of sign in cosz inverts 
the action of the given torque MtA at the DNA level, 
and we see that such a sign alternation actually occurs 
in each interval of length 10.6 bp. 

Conversely, a modification of the linker twist, as can be 
achieved by intercalator enzymes, will modify in a tun- 
able way the geometry of the SH fiber; according to the 
SH chirality, the same change of the linker twist will ei- 
ther condense {AD < 0), either decondense {AD > 0) 
the chromatin fiber. This possible mechanism of conden- 
sation/decondcnsation and its biological relevance will be 
investigated in a following paper (see also Section V.F). 



E. Special geometries 



Ribbon (or zigzag): a first peculiar geometry of the fiber 



is the zigzag or ribbon-like structure, for which the cross- 
linked and roughly toroidal coiling degenerates into a flat 
structure. The number of nucleosomes per turn then 
reaches its minimal value of 2, i.e. 9 = n. In this case, 
represented on Figure 14, 



length change (An^j, 
den configurations. 



1) drives it into sterically forbid- 
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cosz — sin/3sin((f>/2) 
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Nonchiral configurations: the SH chirality changes when 
cos{r]/2) vanishes, i.e. \r]\ = tt; in this case, the link- 
ers cross the SH axis. This configuration is remarkable 
mainly for its mechanical properties. Indeed, the absence 
of chirality amounts, from its very definition, to the van- 
ishing of the twist-stretch coupling: g = 0. Moreover, as 
explained in Section V.D, the sign of cos(?7/2) determines 
whether the linker DNA is unwound (i.e. denatured) or 
on the contrary twisted when the SH is pulled. 

Since sin(7y/2) = 1, the configuration is the limiting 
case opposite to the helical spring, for which sin(?7/2) = 
0. In particular, when the fiber is only pulled {M = 0) 
the twist energy stored in the linker vanishes: all the 
elastic energy of the linker is stored in the bending degree 
of freedom: 
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(22) 



which never occurs in an helical spring. We under- 
line that non-chiral structures (special linker positioning) 
and ribbon-like structure (special nucleosome position- 
ing) should not be confused. 

In fact, this configuration is not rotationally symmet- 
ric, hence has two bend elastic constants. They are com- 
puted within the same lines as above, but without per- 
forming any average over the direction of the bending 
torque Mf,; the identification is done with an asymmetric 
EWLR (4x4 matrix with two bend persistence lengths). 
This yields: 
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Nevertheless, the orientation of the bending torque is not 
controlled at this level, hence the relevant elastic behav- 
ior, to be observed experimentally, is rather that pre- 
dicted after having performed an average over the torque 
directions. This recovers a rotationally symmetric behav- 
ior, characterized by the above given persistence length 
A, as it can be seen directly in the formula: 
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We locate this kind of configurations by determining on 
Figure 9 the linker lengths for which 2-11/9 equals 2. It ap- 
pears (see Figure 10) that cos z then reaches is maximal 
value. Note that a ribbon is chiral, and is not the max- 
imally extended configuration; on the contrary, a small 



Columnar packing: special configurations are obtained 
when ribp is close to an integral multiple of n^^ (for in- 
stance Ubp = 42 bp or 43 bp). In this case, nucleosome 
axes are all parallel hence parallel to the SH axis due 
to the rotational symmetry. We actually check on the 
curves of Figure 10 that (3 — Q for Ubp ~ 42 bp, whatever 
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$. It is noticeable that geometric characteristics as D, z 
and r] are much sensitive to variations of Uhp around this 
value of 42 bp, which reflects in a similar sensitivity of 
the elastic constants, as seen on Figure 11. Fine tuning 
of rifjp and leads to a configuration in which the nu- 
cleosomes are organized in columnar arrays. It has been 
evidenced l^^l that free nucleosomes exhibit a liquid- 
crystal-like nature and tend to form columnar phases, in 
which the nucleosomes are stacked one upon each other 
with their axes parallel to the column. A SH configu- 
ration exhibiting a columnar packing of nucleosomes in 
low-salt conditions is likely to be stabilized by internucle- 
osomal interactions when the ionic force increases, better 
than any other configuration. In presence of linker his- 
tones, the specificity of "columnar" SH is strengthened 
since linker histones too exhibit a propensity to stacking. 

F. Comparison with chromatin fiber pulling 
experiments 

We carry on the discussion by presenting the experi- 
mental validation and prospects of our results. 

In the past few years, a novel experimental method of 
investigating structural and mechanical properties of bio- 
logical macromolecules or complexes, as DNA or chromo- 
somes, has been developped [0 ||ll[] Q |Q. It is based 
on micromanipulations and force measurements on an 
isolated fiber. It belongs to the rapidly expanding field 
of investigations known as "single molecule biophysics" 
p2| . The experiment deviccd in the context of chromatin 
study transposes to a chromatin fiber a methodology first 
implemented with DNA. It consists in pulling the fiber 
under various constraints and in various conditions (salt 
concentration, presence of specific enzymes or chemical 
factors), and to determine the ensuing deformations of 
the fiber. Varying the pulling force yields force-extension 
curves characterizing the elastic response of the fiber. 

Nevertheless, at the chromatin scale, the applied 
stresses are artificial and cannot sensibly refer to an event 
occuring in vivo: contrary to DNA case, an enzyme is not 
large enough to directly handle the 30-nm fiber or to ex- 
perience the chromatin fiber state of strain, for example 
its torsional strain or its curvature (the only in vivo mech- 
anism of that kind, i.e. the action of the mitotic spindle 
on chromosomes occurs at a much larger scale). The 
biochemical processes, for instance binding of a biological 
factor or chemical modification, occur at the elementary 
level (scale of DNA and nucleosomes), whereas the ex- 
periment probes the 30-nm fiber behavior. A mechanical 
modeling relating the DNA scale and the fiber scale is 
thus necessary to exploit all the informations provided 
by single-fiber pulling experiments in terms of biological 
functions. Conversely, these experiments are essential 
to validate the model and the underlying hypotheses on 
structure and interactions, to fit dubious parameters and 
possibly to ask for refinements. The EWLR model allows 



to predict force-extension curves, to be compared with 
those obtained in single-fiber experiments. Our study, 
describing the microscopic root of the EWLR model of 
chromatin, thus provides a bridge between microscopic 
structure and observable properties of the overall fiber. 
In the experimental setup used by Cui and Bustamante 
the fiber, of relaxed length L, is pulled by means of 
optical tweezers. There is no direct contact with the 
fiber ends, which are free to rotate. No torque is ap- 
plied {M = 0) and free twist fiuctuations take place. It 
corresponds to the situation recalled just above. Some 
of force-extension curves that they present are relax- 
ation curves in low salt, for which it is legitimate to ig- 
nore interactions between nucleosomes. Cui and Busta- 
mante fitted these experimental relaxation curves within 
the EWLR model presented in Section IV. A. Recalling 
that their experimental setup allows free twist fluctu- 
ations, the fit involved only two elastic constants: the 
bend persistence length A and an effective stretch mod- 
ulus Jeff — 1 " ksTg^/C. They found A ~ 30 nm and 
7e// ~ 5 pN, values to which we may compare our theo- 
retical predictions. 

The comparison leads to a striking result: we evidence 
on Figure 1 1 that the set of values {A « 30 nm, 7e/ / ~ 5 
pN) is obtained for n^p belonging to a narrow window 
between 42 and 43 bp whatever the value of $ is (we 
checked values of <I> ranging from 130° to 40°). Since n^p 
is then an integral multiple of n^^, this case corresponds 
to special structures of the SH in which the nucleosomes 
pack into columnar arrays (IV. E); only the number of 
columns varies with $ and with the precise modeling of 
the chromatosome, see Figure 3. Pulling the fiber will 
induce a change in $ but not in r, due to the special 
orientation of the nucleosomes; this orientation is thus 
preserved, hence 7 and A will not change: no nonlinear 
effect arises when starting in this special columnar geom- 
etry, hence a fit by an EWLR actually accounts for the 
whole force-extension curve, even in the region where u 
is not small with respect to 1. 

In an experiment performed at higher salt concentra- 
tion (in 40-150 mM NaCl) Cui and Bustamante observed 
a plateau in the force-extension curve at a value Fc be- 
tween 5 pN and 6 pN. They interpret this plateau as 
a structural transition corresponding to the breaking of 
some short-range attractive interactions between chro- 
matosomes, from which follows a dramatic decondensa- 
tion of the fiber, at constant force. Our modeling cannot, 
of course, directly account for this conformational transi- 
tion since internucleosomal interactions are ignored: only 
relaxation curves can be predicted, or curves in low salt 
conditions where interactions are always negligible, even 
at low force, and the chromatin always extended (there 
is besides no plateau in this low-salt case). 

But we claim that our study is specially well-suited 
to investigate the relevant biological question, which is 
rather to determine local mechanisms that could in- 
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duce such a decondensation, i.e. that could create lo- 
cal stresses of large enough strength to break the inter- 
actions. Indeed, a simple criterion of efhciency of bio- 
logical factors (e.g. enzymes) is that the global stress 
resulting from their binding should be larger than the 
critical force Fc measured by Cui and Bustamante. Our 
approach provides a direct method to check this criterion 
since it relates the local strains induced by the binding of 
the enzyme on linker DNA (experimentally measured or 
deduced from a molecular mechanics simulation) first to 
the associated local stresses (equation ^l|); then, by ex- 
tending the computations presented in Section IV, it is 
possible to relate these local stresses to the global stress 
experienced by the chromatin fiber, to be compared to 
the measured value Fc, and to the global deformation 
of the fiber conformation. The implementation of this 
methodology in the case of intercalators, which gives in- 
sights on the in vivo condensation/decondensation mech- 
anisms of the chromatin fiber, will be presented in a sub- 
sequent paper. 

We mention that our work is not the first theoretical, 
model-based approach aiming at accounting for these ex- 
perimental results. 

A first one has been developped by Katritch et al. Q, 
using a Monte Carlo simulation of the fiber. One of the 
interests of such a simulation is to include explicitly ther- 
mal fluctuations but, as they show in their paper, the 
behavior experimentally observed is mainly determinis- 
tic (except at very low force). Moreover, the possibility 
offered by a simulation to treat accurately the interac- 
tions, for instance between nucleosomes, is at the mo- 
ment hampred by the lack of data, preventing to go be- 
yond an effective isotropic model of interaction. We thus 
believe that Monte Carlo simulations are not necessary, 
at least as concerns the problem of reproducing the ex- 
perimentally observed force-extension curves. Indeed, an 
analytic answer is at hand, even with nucleosome interac- 
tions, and with more refined excluded- volume constraints 
(see Figure 7 and the analog that could be constructed 
for the linkers, by unwrapping the cylinder of radius r). 
Simulations would become really useful to handle a more 
detailed description of the chromatosome, for instance 
with an explicit linker histone. 

A second approach has been performed recently by 
Schiessel et al. 0. As discussed below (VII.B), its ba- 
sic step is to relate the strains at the DNA and fiber 
scales (instead of the stresses, as performed here in Sec- 
tion IV). They only managed to determine some elastic 
constants (namely 7 and A) of special geometries, and at 
the expense of some approximations, but they give an- 
alytical expressions directly in terms of the microscopic 
parameters ( here / and $). Moreover, the forces that 
they predict in their theoretical force-extension curves 
are smaller than the results of Cui and Bustamante by a 
noticeable factor of 4. This could be explained by their 
derivation, summing up the bend and twist contributions 
to 7""'^, which is questionable. Indeed, looking at the ba- 



sic example of an helical spring shows that instead of 
writing 7-HAC) = j-\A,C = 0) + -f-^{A = 0,C), 
one should sum up the stretch moduli, according to 
7(A, C) = 7(A, C = 0) + -/{A = 0, C) as in the case 
when springs act in parallel. This leads to a discrepancy 
by a factor of 4 when the "elementary" strech modulus 
j{A, C = 0) and j{A = 0, C) are of the same order. 
Moreover, in the case of the SH, the actual expression - 
see equation (|l^)- shows that such an additive decom- 
position does not exactly hold in the SH case. These 
two objections could explain the discrepancy between our 
value of 5 pN for 7, in agreement with the experimental 
results of Cui and Bustamante, and their value 7 = 1.2 
pN, yet observed for quite similar fiber conformations, 
with crossed linkers (Schiessel et al. considered a struc- 
ture with (in our notations) $ = 85° and r = 36°, which 
corresponds to 42-43 bp in our model). 

VI. BIOLOGICAL DISCUSSION 

A. A tunable, highly sensitive, elastic structure 

We thus evidence a wide range of different elastic be- 
haviors, separated by a minor change in the linker length. 
We suggest that this tunable elasticity might be used as a 
regulatory mechanism during the cell cycle. For instance, 
a slow modulation of the linker lengths, might create dif- 
ferent domains in the chromatin fiber, of much different 
rigidities, and might provide a preliminary underlining of 
transcriptionally active chromatin regions. The response 
of these different regions to a same local stress (protein 
binding, for instance) will be dramatically different. For 
example, intercalation might condense or decondense the 
fiber, according to the sign of dD/dr, i.e. to the chiral- 
ity: within Aubp = 2 bp, opposite consequences will be 
observed. This underlining, inscribed in the very struc- 
ture of the fiber, allows a rapid and selected response to a 
non-specific stress, which might be biologically more rel- 
evant than a mechanism based on enzyme recognition of 
a specific sequence, moreover possibly buried inside the 
fiber. Mechanical sensitivity is likely to provide efficient 
switches for processes occuring at the fiber level. 



The chromatin fiber thus exhibits tunable structure, 
tunable chirality and tunable elastic properties. We sug- 
gest three possible mechanisms to implement the required 
adaptation of the linker length: 

- the first mechanism involves nucleosome displace- 
ment; acetylation of histone tails untightens the DNA 
wrapping around the histone core and presumably allows 
nucleosome mobility. Nevertheless, topological (linking 
number conservation) and mechanical (helical gearing) 
constraints make the motion of the nucleosome far dif- 
ferent from a mere translation along the DNA, and the 
kinematic feasibility deserves to be investigated further; 
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- the second mechanism involves intercalating en- 
zymes, modifying the twist of the linker. A detailed study 
of the interplay between linker intercalation and the chro- 
matin fiber mechanics will be presented in a subsequent 
paper; 

- the best candidate might be a mechanism involving 
the linker histone. Indeed, the value nbp involved here is 
the effective length of the linker, beyond linker histone 
(i.e. outside the chromatosome). This length might be 
tuned by a slight displacement of the linker histone away 
{Auhp < 0) or towards {Auhp > 0) the nucleosome. 

An alternative tuning mechanism lays on the variation of 
<i>, controlled in particular by the presence of linker his- 
tone, salt concentration and histone tails binding affini- 
ties. 



B. A novel chromatin structure 

We evidenced that a value of linker length n^p between 
42 and 43 bp leads to the values A ^ 30 nm and 7 « 5 
pN whatever $ . Since Ubp is then about 4 n°p, the nucle- 
osome axes are all parallel hence parallel to the SH axis 
due to the rotational symmetry. That such organized 
structures lead to ,4 = 30 nm and Jeff = 5 pN, whatever 

is (and even whatever rnuci and H are, as we checked) 
is explained by the fact that the SH elastic properties 
originate from the linker DNA contribution. As shown 
by the computation of SH elastic constants, this contribu- 
tion is mainly fixed by the linker orientation with respect 
to the SH axis (i.e. z). This matching between our pre- 
dictions and experimental results strongly suggests that 
the structure underlying the observed elastic properties is 
a columnar packing of nucleosomes. We thus supplement 
the argument of Yao et al. |^ in favour to rotationally 
phased nucleosomes: the actual relaxed SH structures 
are selected according to their ability to be stabilized 
by internucleosomal interactions when the ionic force in- 
creases. Our claim is supported by the results of Livolant 
et al. 1^ Q, in which it is observed that nucleosomes 
exhibit a liquid-crystal-like nature, leading to a spon- 
taneous columnar ordering. We thus expect native SH 
structures to favour columnar packing of nucleosomes, 
as they correspond to the more robust three-dimensional 
organization of the chromatin fiber. Note that whereas 
a model with straight linkers is compatible with the ob- 
served structure and elasticity of the chromatin fiber at 
low ionic strength (5 mM NaCl), interactions between 
stacked nucleosomes should induce a bending of the link- 
ers at higher ionic strength. More probably, the conflict- 
ing effects consistency between nucleosomes stacking and 
linker stiffness might be reconciled by linker DNA kinks, 
occuring near the entry/exit points and induced by the 
binding oh linker histones HI . 

For $ = 90° (low-salt situation), the degree of com- 
paction 10/Z)(nm) reaches its minimal value for the same 
value rifcp located between 42 and 43 bp. This under- 



lines a key feature of the corresponding configuration: 
at low salt, it is the most extended and rigid configura- 
tion; at the same time, it is the most responsive to salt 
induced compaction. Indeed, as seen on Figure 9, the 
degree of compaction 10/Z3(nm) is minimal at Ubp = 42 
or 43 bp for <I> = 90", and it strongly increases when <I> 
decreases to 50", which is an acknowledged effect of in- 
creasing the salt concentration. Moreover, this configura- 
tion strongly favours a second compaction step, ensured 
by the attractive interactions between nucleosomes (or 
rather chromatosomes) that arise when the nuclesome 
faces are close enough. We in particular recover in this 
scheme the two-stage compaction of the 30-nm chromatin 
fiber observed experimentally [|6| . An insight on this 
interaction-induced compaction can be obtained by set- 
ting the effective parameters r„„ci and H to 0, thus mim- 
icking the enhanced influence of linker histone at high 
salt; in this case, P decreases to about 6 nm, indicating 
that nucleosomes actually stack very closely onto each 
other and lead to a superstable (and presumably rigid) 
fiber. It has been suggested that HI is required not 
so much to get a folded fiber (compact fibers have been 
observed in absence of HI) but to get a properly folded 
fiber. We suggest that HI might be involved in the tuning 
of the effective linker length and twist, actually involved 
in the assembly. In any cases, HI stacking interactions 
favour configurations exhibiting columnar arrays of nu- 
cleosomes. 



VII. CONCLUSION AND PERSPECTIVES 

A. What did we learn from this study? 

The chromatin sensitivity of the fiber response to 
global or local stresses evidenced in this study sheds light 
on the biological interest of the so peculiar and so univer- 
sal assembly of chromatin fiber. It enlights possible rela- 
tions between small-scale structure and gene regulation 
through the fiber mechanical properties. For instance, 
we evidenced that a minor change of ^ or r around some 
"critical" value inverts its chirality. Hence, its response 
to a torsional torque induces either an extension, either a 
contraction of the fiber, and its response to a pulling force 
will either winds or unwinds the SH fiber. Moreover, as 
DNA itself is chiral (right-handed), the response at the 
DNA level is also controlled by the SH chirality: pulling 
the fiber unwinds the DNA if the SH is left-handed. 

The tunable energy partition between twist and bend 
degrees of freedom at the linker DNA level may be of 
biological interest: according to the chromatin config- 
uration, either a twist-sensitive protein will bind onto 
linker DNA, either a protein whose binding is favoured 
by the local curvature of linker DNA will be the adapted 
factor. Hence the mechanical sensitivity can participate 
to biological recognition or specificity. Conversely, we 
read on the associated curves (see Figure 11) whether 
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a twist- modifying protein (intercalator, gyrase), for in- 
stance, may induce a required strain of the SH. 

Comparison with experimental results give clues about 
the much debated chromatin structure [Q. It brings 
about a novel structure: a columnar packing of nucle- 
osomes. We suggest that a possible role of linker hi- 
stone might be to select the proper structure, at low 
salt, by tuning the linker length so as to have nbp/n'^p 
equal to some integer. Also, a tuning of An^p = 1 or 
2 bp might be achieved by enzyme intercalation. When 
ionic strength increases, compaction takes place, first due 
to a decrease of $ (which keeps the nucleosome orienta- 
tions unchanged, along the SH axis); then, when nucleo- 
somes happen to be stacked in columns, the compaction 
is ensured by the interactions between stacked nucleo- 
somes, between stacked linker histones (counterbalancing 
the linker DNA repulsion) or between histone tails and 
nearby linkers or nucleosomes. 

B. A general method for studying elasticity of linear 
complex fibers 

The problem of relating the elastic coefficients of the 
chromatin fiber to the geometric and elastic properties 
of the underlying "microscopic" structure (assembly of 
nucleosomes and linkers) is reminiscent of similar works 
performed for DNA by Marko and Siggia and O'Hern 
et al. [|2|. They described the dsDNA at two levels: 
as an helical coiling at small scale and as an EWLR at 
a slightly larger scale. They derived similar formulas re- 
lating the elastic coefficients (C, A, jdna and goNA) to 
the geometric and elastic parameters of the underlying 
helical model. 

We point that their computation rests on the rela- 
tion between the "microscopic" strains and the dsDNA 
strains, plugged into the equality of the EWLR free en- 
ergy and the free energy computed within the micro- 
scopic model. Generally, numerous sets of microscopic 
strains achieve the macroscopic strains, but only the set 
of lowest free energy yields the actual free energy of the 
EWLR. Relating properly microscopic and macroscopic 
strains thus requires to minimize the small-scale free en- 
ergy, given the macroscopic strains. Performing this min- 
imzation is nothing but writing the conditions for the 
local equilibrium of the assembly. 

When the microscopic model is homogeneous, the con- 
ditions for local equilibrium simply expresses in the ex- 
tensivity of the strains. The uniformity of the local strain 
densities thus allows to relate them to global strains with- 
out an explicit minimization (think to identical springs 
in series). This works for DNA but not in the case 
of chromatin. The discrete and complex nature of the 
chromatin assembly leads to difficult and cumbersome 
computations in order to determine the linker shape, as 
it can be seen in the work of Schiessel et al. Q]. Moreover, 
this approach fails to give an analytical solution except 



when the relaxed fiber exhibits a special geometry, for 
example a ribbon-like flat structure. 

Determining the conditions for local equilibrium is pre- 
cisely what is done, more straighforwardly, in our ap- 
proach. We indeed write equilibrium equations given the 
global stresses and solve them to get the local stresses 
arising in each point of the assembly (at equilibrium un- 
der the given global stresses). Relating the stresses at 
the microscopic and at the fiber levels follows from the 
basic principles of classical mechanics; this method ap- 
pears to be at the same time more simple and more 
easily generalized. It is in fact the only way to bridge 
the linker elasticity to the SH elasticity in any geometry. 
Moreover, it extends to more complex situations as inter- 
calated linker DNA or more generally, situations where 
forces and torques are applied at the DNA level. 

We thus underline that the proper method to express 
the elastic properties of an assembly as a function of the 
elastic properties of the basic elements is to relate the 
global stresses, applied to the assembly, and the local 
stresses experienced by the lower scale elements. This 
relation can be used in both ways, to investigate 

- either small-scale repercussions of a global stress, 
hence how a global stress (as those applied in microma- 
nipulations) can be used to probe the fiber at the ele- 
mentary level, 

~ either large-scale response to local stresses, hence 
how biological factors binding on the linker DNA could 
induce major structural and conformational changes in 
the overall fiber. 



C. Biological perspectives 

Our study underlines that the mechanical properties 
of special structures, selected according to the phasing of 
the nucleosomes therein, might be involved as a regula- 
tory factor in the chromatin biological function. 

Having modelled the 30-nm fiber as a EWLR, with ex- 
plicit values of the elastic coefficients, a natural extension 
of our study is to consider higher levels of organization: 
a plectonemic coiling, leading to a 60-nm fiber, or an he- 
lical coiling, whose elastic properties follow from classical 
spring mechanics. The question is then to unravel the im- 
plications of the chromatin structure and its elastic prop- 
erties on the higher levels of organization. A mechanical 
approach similar to that implemented in this paper is es- 
sential to bridge electron-microscopy structural observa- 
tions evidencing fiber-like objects at higher levels (60-nm 
fiber, "chromonema" fiber of diameter 100-130 nm ]5^ ] 
and chromosome) and experimental results on chromo- 
some elasticity obtained by pulling a single chromosome 
p3[ iQ. The challenge is to understand the mechanics 
of the chromosome and its involvment in the biological 
functioning of the chromosome throughout the cell cycle. 

Another direction in which to exploit the results of the 
present paper is to determine the stresses that can be 
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exerted at the fiber scale by groove-binding proteins or 
intercalators when they bind onto hnker DNA. For in- 
stance investigate whether a local decondensation of the 
fiber might be induced by intercalating enzymes and con- 
trolled by linker lengths. More generally, our approach 
is a privileged tool to investigate the action at the fiber 
level of small-scale biochemical stresses (protein binding, 
histone tail acetylation) , then to describe how they can 
act as mechanical switches and exploit the tunable elas- 
ticity of the fiber into regulatory schemes. 
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VIII. CAPTIONS 



Figure 1: Geometric description of linker DNA (section 
II. A). Above: the arrows indicate the orientation of the 
two strands, dehmiting the minor groove (fiUed in grey) 
and the major groove. The vector t(s) is tangent to the 
minor groove and rotates around the local axis u (straight 
in a relaxed linker) with an angular rate oj^ — 2Tr/l^; 
the two vectors make a constant angle 90° — aoNA ~ 
62". Below: tranverse view; the bold contour locates the 
minor groove. 

Figure 2: Nucleosome modeling and assembly of the 
chromatin fiber (section II. B). Left above: view in per- 
spective of the nucleosome j; a is the slope of the 
DNA left-handed wrapping around the histone core {a ~ 
4.47°). The dyad axis Dj is orthogonal to the nucleosome 
axis Nj. Left below: view from above (projection on a 
plane orthogonal to the nucleosome axis. Right: vectors 
Nj, Uj and tj{l) in the plane tangent to the nucleosome 
at Ej. 

Figure 3a: Various SH structures obtained when the pa- 
rameter I (equivalently 7if,p) is varied, for $ = 90°. Note 
the columnar packing (/3 « 0) obtained for ribp between 
42 bp and 43 bp, and the ribbon-like structure (two nu- 
cleosomes per turn, 6 = tt) obtained for Ubp = 47 bp, 
whatever the value of In between , the fiber exhibits 
cross-linked configurations. In each period of length 10.6 
bp, an interval (depending of $) of n^p-values is forbid- 
den as it corresponds to self-overlapping configurations 
(here around 7if,p = 38 bp) 

Figure 3b: Various SH structures obtained when the 
parameter / (equivalently Ubp) is varied, for $ = 50°. 
Configurations around Ubp = 45 bp are forbidden (steric 
hindrance) . 

Figure 4: Left: SH from above (projection in a plane 
orthogonal to the SH axis A, here for (3 — 90°. We 
define R as the distance from any nucleosome center Gj 
to the SH axis (measured along the dyad axis Dj), r 
as the distance of any entry (or exit) point to the SH 
axis A; the "excluded-volume" radius Rsh of the SH, i.e. 
the radius of the cylinder of axis A containing the whole 
fiber (including the nucleosomes) satisfies Rsh = R + 
fnuci due to the peculiar orientation of the nucleosomes 
(the axes Nj are tangent to the cylinder of axis A and 
radius R whatever 0). We also introduce the angle 6 
of the rotation around the SH axis A, transforming the 
projection of a nucleosome into the following one, and the 
angle rj between the projections of the vectors relating 
Sj-i and Ej to the SH axis. This angle rj is chosen in 
[0, 27r[. It satisfies Zsinz — 2rsin(?7/2). Right: projection 
of the first nucleosome in the plane orthogonal to its axis 



Ni. Note that the nucleosome dyad axis is orthogonal 
both to the nucleosome axis (nucleosome symmetry) and 
to the SH axis (symmetry of the assembly). 

Figure 5: Left: Location of the nucleosome axis Nj in 
the frame spanned by the SH axis A and the nucleosome 
j dyad axis Dj ; it lays in the plane orthogonal to Dj and 
makes an angle (3 G [0, 27r] with respect to A. This angle, 
varying with <I> and I, satisfies cos/3 = Nj.A. Right: 
Location of the linker Uj in the frame spanned by the 
SH axis A and the nucleosome dyad axis Dj] the linker 
makes a constant angle ^ with the dyad axis. 

Figure 6: Front view of the SH projected in a plane 
spanned by the SH axis A and the linker axis Uj. The 
SH axis is actually ahead of the linker, which lays in the 
figure plane; Ej^i, Gj-i, Sj and Gj are projections on 
the figure plane. Only the track of the nucleosomes is in- 
dicated: neither their axes nor their dyad axes belong to 
the drawing plane. D is the distance between two succes- 
sive nucleosomes along the SH axis A, whose orientation 
is chosen so that D > 0. We denote z the angle between 
the linker and the SH axis; it satisfies Uj.A = cosz and 
by convention z S [0,7r[ hence sinz > 0. 

Figure 7: Unwrapping of the cylinder of radius R (con- 
taining the nucleosome centers). The construction only 
requires to know the values of f3, 6, D and the nucleosome 
dimensions rnuci and H. This drawing can be conve- 
niently used to check the excluded-volume constraint be- 
tween any two nucleosomes. Linkers, crossing the plane 
of the drawing, are not represented, nor the entry and 
exit point (located outside the drawing plane). 

Figure 8: Fiber actual diameter 2Rsh {dashed upper 
line, curve (1)) and pitch P [solid line, curve (2)) in 
nm, versus linker length (in bp), for $ = 90° {left) and 
<& = 50° [right). As what matters is the twist angle 
T — 2Trnbp/n'^p, the number Ubp of bp is allowed to vary 
continuously. The thin line (3) represents the distance 
D between two nucleosomes along the SH axis. The line 
made of crosses (4) represents the lower bound Pc on 
the pitch (equation ||) . The configurations with P > Pc 
satisfy excluded-volume constraints. 

Figure 9: [Solid line) compaction rate 10/-D(nm) (num- 
ber of nucleosomes per 10 nm of fiber) and [crosses) num- 
ber 2tt/0 of nucleosomes per turn versus linker length (in 
bp), for $ = 90° [left) and $ = 50° {right). Two nucleo- 
somes per turn correspond to ribbon-like configurations, 
whereas the maximal number is reached for columnar 
packing {ribp = 43 bp). 

Figure 10: Cos(?7/2) [solid line, (1)), cos/3 [dashed line 
(2)) and cosz [crosses (3)) versus linker length (in bp). 
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for $ = 90° (left) and $ = 50° (right). Cos{r]/2) = cor- 
responds to a change in the SH chirahty; cos /3 = 1 cor- 
responds to columnar packing in which the nucleosomes 
axes are parallel to the SH axis: cosz is then small and 
the linkers are almost orthogonal to the axis (the differ- 
ence is a « 4.5° where tana is the slope of the linkers on 
the nucleosomes). 

Figure 11: Above: Twist persistence length C {dashed 
line, curve (IJ), bend persistence length A {solid line, 
curve (2)) and SH pitch P {lower thick line, curve (3)) 
in nm versus linker length (in bp, according to I — 
fubp/nlp) for $ = 90° (left) and $ = 50° (right). Note 
that they all vanish at a critical value ntp^c, strongly de- 
pending on $ and on the detailed modeling of the nucle- 
osome (or chromatosome) . The dotted line {curve (4)) 
represents the lower bound Pc (equation |^) indicating 
the values of nbp forbidden by steric hindrance. Middle 
stretch modulus 7 {solid line, curve (IJ) and effective 
stretch modulus 7e// — 1 — ksTg^ /C {dashed line, curve 
(2)) in pN for $ = 90° (left) and $ = 50° (right). Be- 
low: twist-stretch coupling g (no dimension) versus linker 
length (in bp) for $ = 90° (left) and $ = 50° (right). 

Figure 12: these figures represent the relative contri- 
bution (in %) of bending energy {curves (IJ) and twist 
energy {curves (2)) to the total elastic energy stored in 
a linker (within the chromatin superhelix), versus linker 
length (in bp), in various instances: {above) when a force 
F = FA is applied along the SH axisA but without any 
torque, {middle) when a torsional torque Mt — MtA is 
applied but no force, {below) when a flexural torque ( 
Mb. A — 0) is applied but no force. Two values of the an- 
gle $ have been investigated: <i> — 90° {left) or $ = 50° 
{right). The almost flat line, below 2 %, in the two upper 
figures, corresponds to the fraction of energy correspond- 
ing to stretch deformations; it is obviously negligible. 

We compare these curves with those obtained for an 
helical coiling of DNA (section V.B) having the same val- 
ues of z and r {crosses, curve (3) for the twist energy and 
diamonds, curve (4-) for the bending energy). When ap- 
plying only a pulling force, the partitions strongly differ; 
by contrast, the partitions observed when only a torque 
(either Mt, either Mb) is applied are exactly identical 
(the curves superimpose). 

Figure 13: Twist persistence length C {dashed line, 
curve (1)), bend persistence length A {solid line, curve 
(2)) in nm of an helical spring versus linker length (in 
bp) with the values of z and r equal to those of the SH, 
respectively for $ = 90° (left) and $ = 50° (right). This 
shows that the decrease of the persistence lengths and 
their vanishing in special values ribp^c is a geometric ef- 
fect: it is due to an accumulation of turns on a small 
distance along the axis, allowing a dramatic unfolding. 



Figure 14: Ribbon- like structure of the SH {9 = n), here 
projected in the plane {D,N A D); indeed, when 9 = tt, 
all the planes {Dj,Nj ADj) coincide. The SH axis makes 
an angle {n/2 — 13) quite small with this plane (7r/2 — = 
a K, 4.47°) the nucleosomes are almost orthogonal to A. 
Obviously, Dsui(3 = /sin($/2). 
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